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1 Introduction 



The one-dimensional asymmetric exclusion process (ASEP) is a lattice model which describes par- 
ticles hopping in a preferred direction with stochastic dynamics and hard core exclusion. It has 
been used to describe various phenomena like growth processes Jl], polymerization |J, and traf- 
fic flow UH . Although the ASEP is one of the simplest stochastic many-particle models it shows a 
rich behavior ranging from phase transitions [||, ||| to the formation of shocks M. During the last 
few years exact solutions have been found for the stationary under various boundary conditions. 
In particular the ASEP with random sequential dynamics and external particle input and output 
at the boundaries has been analyzed in detail || |8[ ||, |i~0f. There are only few exact results for 
dynamic properties (relaxational behavior) of the ASEP flllfl . 

A mathematical tool which has proved to be very useful for the study of stochastic lattice 
models is the so-called matrix-product state (MPS) technique (l^]-[Q,[||]-[|l^]. Matrix-product 
states can be understood as a generalization of ordinary factorizable states with a product measure 
in which numbers are replaced by non-commutative object. This allows the original problem to 
be reformulated in terms of algebraic relations. The main advantage of this technique is that 
once a nontrivial representation of the algebra is available, physical quantities like density-density 
correlation functions can be computed very easily. In 1993 Derrida et. al. showed that the stationary 
state of the ASEP with external particle input and output can be written as a MPS. Since then 



matrix-product methods have been applied to various other problems, e.g. to three-state models [14 



15], excited states in the ASEP [16], certain reaction-diffusion models [17] and models with quasi- 
parallel updates |l8|]-[|(|. Although a lot of work on MPS's has been done, the full range of 
applications is not yet knownPj. 

So far the matrix-product technique has been applied mostly to systems which are homogeneous 
in the bulk. It is therefore interesting to investigate whether this method can also be applied to 
systems with impurities. In case of the ASEP two kinds of impurities are discussed in the literature. 
The first kind of impurities are defective particles which jump with a rate less than that of other 



particles [15, 22, E3]. Such 'moving' impurities can be visualized as slow cars on a motor way 



which (for sufficiently high particle density) induce traffic jams - a phenomenon which has been 
related recently to Bose-Einstein condensation [23|. To this problem the matrix-product technique 
has been applied successfully [15, 22 1. The second kind of impurities are stationary defects which 



can be introduced by lowering the hopping rate at specific bonds. As in the previous case, such 



impurities can provoke the formation of shocks [24]. Defects of this kind play an important role 



in traffic flow applications wherefore exact solutions are particularly interesting. However, for 
some reason ASEP's with stationary defects are much harder to solve than systems with moving 
impurities. Exact results exist only in case of an exclusion model with parallel dynamics and 
deterministic hopping in the bulk [25]-[p(|. The case of stochastic hopping or random sequential 



dynamics has not yet been solved exactly. Moreover matrix-product techniques have not been 
applied to this type of problems. 

In this paper we consider the stationary state of the ASEP on a ring with various kinds of 
parallel dynamics, deterministic hopping in the bulk and a single stationary defect. Using Bethe- 



Ansatz techniques, this problem was first solved in Ref. [26] in case of alternating parallel updates 



3 As it has been shown recently, for every one-dimensional stochastic lattice model with local random sequential 
dynamics the eigenstates of the Hamiltonian can be written as MPS's [^lj. However, matrix representations which 
are useful for practical purposes seem to be limited to a few models. 
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Figure 1: The ASEP on a ring with L = 2N = 8 sites. In the bulk particles hop deterministi- 
cally in clockwise direction. A defect is introduced between sites L and 1 where particles hop with 
probability q. 

on two sublattices. The solution, however, is fairly complicated and amounts to a list of rules for 
an operative construction of the stationary state. In the present work we solve the problem in a 
more compact way by using the matrix-product technique. Moreover this formalism allows models 
with different update sequences to be solved within the same framework. In order to demonstrate 
this, we consider two different update sequences . Our main intention, however, is to outline a 
method which might help to solve the problem of stochastic hopping in presence of a defect. 

The paper is organized as follows. In Sect. [2] we define the ASEP on a ring with a defect and 
introduce two different dynamical rules. In Sect. ^ we show that the two-state model on a ring 
can be reformulated as a four-state model on a linear chain. Using this mapping we formulate a 
matrix-product ansatz (see Sect. leading to a set of algebraic equations which turns out to be 
the same for both update sequences. In Sect. [| we give simple two-dimensional representations 
of this algebra. Using these representations, some physical quantities are derived and discussed in 
Sect. ^ Finally in Sect. [?] we summarize our results and discuss possible generalizations. 

2 The model 

The exclusion model we are going to study is defined as follows: Particles of one species move in 
clockwise direction on a one-dimensional ring with an even number of sites L = 2N. Each site 
can be either free or occupied by one particle. The time evolution is discrete and defined by a 
certain sequence of parallel updates to be specified below. In the bulk of the chain these updates 
are deterministic, i.e. in each time step the particles move forward provided that the following site 
is empty. Between two particular sites (by convention sites L and 1) a defective bond is inserted 
where the particles hop stochastically with a given probability < q < 1 (see Fig. 1). 

In what follows we will consider two different update sequences. The first one we call sym- 
metrized sequential dynamics: In each time step the first update takes place at the sites (N,N + 1) 
just opposite the defect. Then a sequence of pairwise updates follows. These updates are ar- 
ranged symmetrically with respect to the defect. The first pair to be updated is (N — 1,N) and 
(N + 1,N + 2), followed by (N — 2,N — 1) and (N + 2, N + 3) up to (1, 2) and (L - 1, L). The last 
pair of sites to be updated is the defect (L, 1). As we will see below, this is the simplest dynamical 



rule in terms of the matrix formalism. 

The other dynamical rule we will use is called sublattice-parallel dynamics. Here we have to 
assume that N = L/2 is an even number. Each time step consists of two separate half time steps. 
In the first half time step particles located at odd sites move one step in clockwise direction pro- 
vided that the next site is empty. Then in the second half time step particles at even sites move 
forward in the same way except for the defect where such moves take place with probability q. This 
dynamical rule was introduced in Ref. for solving the deterministic ASEP with a defect. A 



similar dynamics was also used in Refs. [25, 18] in case of the ASEP with external particle input 
and output. 

Symmetrized sequential and sublattice-parallel dynamics can be cast in a more formal way as 
follows. Let Tj = 0, 1 be the occupation number at site j and consider the space of all configurations 
in a canonical configuration basis. Let T (T^) be the two-site hopping matrix in the bulk (at the 
defect): 
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Then the transfer matrix T seq for symmetrized sequential dynamics reads 



seq 



T L \ (Tl^TL-I.l) (^2,3^L-2,L-l) • • • (Tn-1,nTt 



l,N 2 N+l,N+2) Tn,N+1 , 



where the indices indicate pairs of sites on which T acts. The transfer matrix for sublattice-parallel 
updates factorizes into two factors 

(3) 



rp _ _ rp(2) rp(l) 

-Lpar — - L par x par 



which are given by 
T« = 

^par 



T (2) 
^par 



{Tl,2^L-l,L) (^3,4^1-3,1,-2) • • • (Xn-1,nTn+1,N+2) , 



(4) 



N+2,N+3)?N,N+1 



Notice that both transfer matrices T par and T seq differ only in the sequence of their (noncommu- 
tative) factors. There are, however, many other transfer matrices with different update sequences. 
As we will see below, the matrix-product formalism allows to compute ground states for various 
dynamical rules from the same matrix representation just by rearranging the factors in the ma- 
trix product. This is an important advantage of the matrix-product method compared to direct 
techniques as presented in Ref. [pH|. 



3 Reformulation as a four state model 



Although the defect breaks translational invariance, the model is still symmetric under reflections 
with respect to the origin combined with particle-hole symmetry. Since we expect correlations in 
the stationary state to be subject to this symmetry, it is natural to introduce a matrix ansatz which 
exploits this symmetry. This can be done by grouping pairs of sites together which are located 
symmetrically with respect to the defect. This mapping defines a four-state model on a linear 
chain with iV sites and closed boundaries which is equivalent to the original one (see Fig. ^). In 
this language the two-site hopping matrix T^i m * ne original notation is equivalent to an one-site 
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Figure 2: By grouping pairs of sites, the two-state ASEP on a ring can be regarded as a four-state 
model on a linear chain. 



operator C\ at the left boundary acting in a four-dimensional space. Similarly 7zv,jv+1 corresponds 
to an one-site operator TZn at the right boundary. The block-spin variables U{ = (rz-i+i , t*) 
can take four different values {0, 1, 2, 3} = {(0, 0), (0, 1), (1, 0), (1, 1)}. In this basis the boundary 
hopping matrices read 
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(5) 



In the bulk pairs of opposite hopping matrices TL-i,L-i+l) are grouped together resulting in 

a two-site operator «Sj j-fi which is given by a 16 x 16 matrix: 
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In terms of these operators, the transfer matrix for sequential dynamics (||) takes the simple form 



- seq 



C-l <Sl,2 ^2,3 • • • Sn-1,N T^n ■ (7) 
Similarly the transfer matrix for symmetrized sublattice-parallel dynamics (|j) can be written as 

^par = <Sl,2 ^3,4 • • • Sn-1,N j ^par = A ^2,3 ^4,5 • • • <?7V-2,iV-l T^N ■ (8) 

So far, we have only rewritten the original model on a ring as a four state model on a linear chain. 
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4 The matrix-product formalism 



As we will show, the four-state formulation allows to write the stationary state of the asymmetric 
exclusion model with a defect as a matrix-product state. Let us define noncommutative operators 
(matrices) Aq, A±, A2, A3 and Bq, B\, B2, B3 as well as bra and ket vectors (W\ and |V) in some 
auxiliary space and demand that they obey the relations 

S(A®B) = (B®A) (9) 
KA\V) = B\V) (10) 
(W\CB = (W\A . (11) 

where A = (Aq, A±, A2, A3), B = (Bq, B\,B<i, B3) and '(g)' denotes the tensor product in configura- 
tion space. This algebra can be used to construct the ground state of the stochastic model we are 
studying. Let us first consider the case of symmetrized sequential dynamics. As we will show below, 
the probability P seq (a\, . . . , ajv) to find the stationary system in the configuration {a±, . . . , a^} can 
be written as 

P seq (a 1 ,...,a N ) = -= m (W\A tT1 A tT2 ...A tTN \V), (12) 

N 

where Z^ q is a normalization constant which is given by 

Z s N eq = (W\C N \V) with C = A + A 1 + A 2 + A 3 . (13) 

For example, the probability to find the four-state chain with N = 4 sites in the configuration 
{1, 0, 3, 2} = {(0, 1), (0, 0), (1, 1), (1,0)} is given by P seq (l, 0, 3, 2) = (W\A 1 A A 3 A 2 \V) / \W\C^\V) . 

The mechanism ensuring the stationarity of these probabilities is precisely the one introduced in 
Ref. |9| and works as follows: Formally Eq. ([[D may be written as \P seq ) = (ZP)' 1 (W\A® N \V). 
Applying the transfer matrix T seq = C\ S\ : 2 52,3 • • • Sn-i,n T^n to this state the matrix 7Z at 
the right end of the chain generates the vector B = (Bq, B\, B2, B3) at the iV-th position in the 
product (cf. Eq. (0)). Then, by successively applying S, the generated vector B is commuted 
to the left (cf. Eq. (g)). Finally, when reaching the left boundary, the vector B is turned into A 
by the action of the defect hopping matrix C (cf. Eq. (pi])). Consequently, the application of the 
transfer matrix T seq on the state \P se q) results in the same state |-P se <j) • Thus the state \P se q) is a 
stationary state of the transfer matrix ([?]), i.e. of the transfer matrix (|2|). 

For sublattice-parallel dynamics we have to use a slightly different ansatz which involves alter- 
nating matrices. We write the stationary state as 

P par (a 1 ,...,a N ) = 1 ^ar(W\(A ai B (72 )(A (73 B a4 )...(A aN _ 1 B aN )\V), (14) 



7 P 
N 



where 



Z p N ar = {W\(Y,A l Y,B j ) N/2 \V). (15) 

i=i j=i 

For this type of dynamics the mechanism ensuring the stationarity of the probabilities (|l~4| ) is 
precisely the one proposed in Ref. Jl8| . Again we use use the algebraic relations (^)-(ll). Writing 
\Ppar) = (-^jv"*) -1 (W\ (A (g> B)® N I 2 \V) one can easily verify that in each half time step the vectors 
A and B are exchanged: 

T£l[(A®B)® N / 2 ] = (B®A) m / 2 , (16) 
(W\T^H[(B <g) A)® N I 2 ]\V) = {W\(A® B)® N/2 \V) . (17) 



5 



Therefore the state \P pa r) is invariant under application of the transfer matrix T par . 

It should be emphasized that both symmetrized sequential and sublattice-parallel dynamics 
use the same algebra (^)-(ll). One can also consider various other update sequences leading to 



(H)"(0)- For example, if updating is started in the middle of the system and proceeds to both 
ends from there, the stationary state can be written as (W\ A® N B® N \ V) with the Ai , Bi and 
(W\ , \V) obeying In other words, having determined a nontrivial representation of the 

algebra given above, one can immediately compute stationary states for various update sequences. 



5 Representation of the algebra 



Eqs. (P)-(|lT|) define an algebra of eight objects Ao, A\, A2, A 3 , Bq, B\, B2, B 3 . There are sixteen 
bulk equations 

[Ai,Bi] = (i = 0...3) 

B A 2 = B A 3 = B X A Q = B X A 2 = B X A 3 = B 3 A = B 3 A 2 = (18) 

AtBj+AjBi = BiAj e {(0,1), (2,0), (2,3), (3,1)} 

^0-63 + AiB 2 + A 2 Bi + A 3 B = B 2 A X , 

four equations at the defect (left boundary of the four-state model) 

(W\B = (W\A 

(W\(B l + qB 2 ) = {W\A X (19) 

{W\{l-q)B 2 = {W\A 2 

(W\B 3 = (W\A 3 , 

and four equations opposite of the defect (right boundary): 

A \V) = Bq\V) 

= B^V) (20) 

(A 1 + A 2 )\V) = B 2 \V) 

A 3 \V) = B 3 \V). 



Notice that the algebra is invariant under the replacement J27 

^lo^A^o, Bo^XB , A 3 — > A -1 A 3 , B 3 — > A -1 B 3 , (21) 
where A is some number. 



The algebra ([Lq)-(|20|) has a complex structure and thus it is practically impossible to determine 
its representations directly. In order to simplify the problem, we therefore use the following ansatz: 

B = A , B 1 = A 1 -1, B 2 = A 2 + 1, B 3 = A 3 . (22) 



A similar ansatz has been used in Refs. [19, ^] in order to relate sublattice-parallel and random- 



sequential updates in the ASEP with external particle input and output. It is obvious that Eq. (p2j) 
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imposes strong constraints and therewith reduces the space of solutions. However, it turns out that 
it still includes nontrivial physical representations. 

Inserting Eq. fl2^) into Eqs. (|l|)-(|20|), we obtain a reduced algebra of four objects. It consists 
of seven bulk equations 

A X A = A , A X A 3 = A 3 , A X A 2 = A 2 , A A 2 = A 3 A 2 = A A 3 = A 3 A = (23) 
and two boundary conditions: 

(W\A 2 =^—^-(W\, A l \V) = \V). (24) 

This algebra is much simpler and can be analyzed systematically on a computer. In fact, we found 
the following two-dimensional representation: 











\ 




i ■ 




h- 1 



A) = A n 1 . M= \ n i I, A 2 = 1 n n . ^3 = 0, 




W = (l,l), 1*0 = ft) > ( 2 5) 



where A is a free parameter. Algebraically A is related to the invariance of the algebra in Eq. 
Physically it is related to the conservation of the number of particles. In fact, as will be discussed 
in the next section, the ground states ( |i~2| ) and ( |l~4|) describe grand-canonical ensembles of systems 
with different particle numbers where the parameter A plays the role of a fugacity 



As can be verified easily, the canonical ensemble described by the above representation includes 
all sectors with M < N = L/2 particles (if there were more than L/2 particles, at least one of 
the block spins ui = (7x_j+i , n) would be in the state (1,1). Because of A 3 = this implies that 
the corresponding matrix product vanishes). However, the algebra (p3j)-(|24|) is invariant under the 
exchange of the matrices Aq «-> A 3 which immediately gives a second set of representations for 
systems with more than L/2 particles. Physically this invariance is related to the particle-hole 
symmetry in the ASEP. 



6 Some physical quantities for the model with symmetrized se- 
quential dynamics 

6.1 Micro-canonical ensemble: Density-profile for a fixed particle number 

In order to derive physical quantities for a fixed particle-number, we have to project the stationary 
state onto a sector with a specified number of particles M. A projection formalism for matrix- 
product states has been introduced recently in Ref. J22| |. We are now going to apply this formalism 
to the present model with sequential dynamics (the case of sublattice-parallel dynamics can be 
treated similarly). Because of particle-hole symmetry, we restrict ourselves to less than half-filled 
systems where M < N = L/2. 

Let x( m ; • • • > a n) be a function which is 1 if the sites 1 . . . n are occupied by m particles and 
otherwise: 

n 

x(.m;a 1 ,...,a n ) = 8{n - m + ^(<W - ^,0)) • (26) 
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Then the probability to find the model with sequential updates and M particles in the configuration 
{at, ... ,a N } is given by 

P^{a u ....a N ) = — ^- x{M;a u ...,<r N ) (W\ A ai . . . A aN \V) , (27) 

N,M 

where Z^ q M is the normalization constant restricted to the M-particle sector: 

Z N q M = (W\G NM \V) (28) 
G N , M = x(M;a 1 ,...,a N )A ai ...A aN (29) 

ffi, ...,er/v 

The expression GV,m is the sum of all products of N matrices with M particles. By definition, two 
of these objects can be combined by convolution: 

M 

Gn!+N 2 ,M = GWi,j Gn 2 ,m-j ■ (30) 

3=0 

Physical quantities can be expressed in terms of combinations of these objects. For example the 
density c{x) to find a particle at site x in the original model on the ring with L = 2N sites and 
M < N particles is given by 



1 f Efif" 1 '^ (W\G x - ltj A x G n ^m- 3 -i \V) if x < L/2 
7j "M \ Y™^ L ~ xM ~ l) (W\Gl- x ,3 A 2 G^jv-i.M-i-i \V) if x > L/2 



Similar expressions exist for the current and higher correlation functions. Once G n ^ m is known, all 
these quantities can be computed immediately. 

The expression G n ^ m can be defined recursively by G n>m = if m < or m > n, Gqq = 1, and 

G n ,m — ^0 G n — l )J7 j 

+ (Ai + A 2 ) G n _i jm _i + A 3 G n -\.m-2 ■ (32) 

However, instead of solving this recurrence relation algebraically, it is much simpler to use directly 
the representation fl25|) . We obtain 

where 

/ \ m / \ 

n \ „ \ / n 



Using this result we obtain the normalization ( |2"8| ) 

q j q M -i = 2q- M (1 - q) M ~ N h_ q (N - M, M + 1) , (35) 



M / 

Z S N q M = 2<T M E 

3=0 V 



N 
M-j 



where I z (n,m) is the regularized incomplete beta function. Now the particle density profile ([H]) 
can be computed easily. Using that A\ = 1 and (W\G^'L i A2 = (W\5 ntJn (l — q)q~ n ~ l we get 

c(x) = -J- J Z " 1 ' M - 1 ^ X " (36) 

( } Z Zf 1(1 - q) <f~^ C™, if - > 1 ^ 
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This formula holds for less than half filling M < N = L/2 (the case of more than L/2 particles 
is related by particle- hole symmetry). Writing x = L — y + 1 and keeping the particle density 
p = M/N <C q fixed one can now derive the asymptotic expression for the particle density in front 
of the defect in the thermodynamic limit: 



Jim c{L-y+l) = (1 - q) (37) 
N— >oo Vfir/ 



A similar result was derived in Ref. [26] for sublattice-parallel updating. 



6.2 Grand-canonical ensemble: Correlations for a large lattice 

Let us consider a grand-canonical ensemble of systems with L S> 1 and an average particle number 
pL where p is some given density. Its probability distribution is given by Eq. (jl^) with the matri- 
ces (p5|). The 'fugacity'-parameter A has to be chosen such that the particle-density is equal to p , 
i.e. it has to solve the following equation: 

1 N 

2Np = -agj £ (W\ C*- 1 (A, +A 2 + 2A 3 ) C N ^ \V) . (38) 

Z N i= i 

Like in the previous section we restrict to particle numbers less than N = L/2 (Results for systems 
with higher particle numbers are obtained exploiting the particle-hole symmetry). Because of 
C = A Q + A\ + A 2 + A 3 and ^3 = Eq. (|||) can be written as 

1 N 

1 - 2p = -—sat F N,i (39) 

iV Z N i= i 

with F Nii = (WlC^AoC^lV) . 

The computation of the matrix-elements Z^ q and Fjf i is done most easily in a representation where 
C is diagonal. Such a representation can be obtained by means of a similarity transformation from 
the representation (|25|). We will use the following matrices and vectors 

c-K 1 ,!,), ^ V <«) 



to compute Ztf 9 and Fjy j as 



ryseq 
N 



q + Xq 



1 + A y u ^0 1 



- - { Xq (1 + A) N + (? - 1) } (41) 



The large- iV asymptotics of these expressions depend on the magnitudes of the terms 1 + A and q^ 1 . 
Let us consider the case 1 + A > q- 1 with l + X-q- 1 = 0{N°) first. Inserting Eqs. Q) and (@) 
into Eq. (|39|) and approximating for iV 3> 1 results in 

A = ^. (43) 
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This equation relates the parameter A to the density p under the condition 1 + A > q 1 , i.e. for 

P < | • (44) 



Eq. (43) completes the description of the system for the case p < q/2 . One can now easily compute 
the density-profile or correlation functions. Noticing that C plays the role of a transfer matrix one 
reads from Eq. (|io| ) that all correlations decay exponentially (see e.g. Ref. |fT?| l) with a correlation 
length £ = {log(l + A)g} _1 . Using the explicit expression (|43|) for A we find 



log^-j 1 . (iV»l, p<q/2) (45) 

This result is the same as the one obtained in Ref. [26] for the model with sublattice-parallel 
dynamics. 

In the case p > q/2 the above assumptions on A do not lead to any result. In order to obtain a 
finite density we have to assume A to be of the form 

A = -1 + q- 1 ( 1 + £) (46) 

with some number A' which is independent on iV . Computing again the large-iV asymptotics of 
Eqs. ([II]) and ( f42| ) and inserting it into Eq. (|39| ) results in the transcendental equation 

_ 1 + (A' - l)e V 
1-q A'(e A '-l) [ ' 



which is soluble for densities p > q/2. Eqs. ( f46|) and fl47|) relate the parameter A to the density p . 
The correlation length, which is given by £ = {log(l + X)q}~ 1 ( £ = — {log{\ + \)q\~ l ) for A' > 
( A' < ) takes now the form 

£ ~ 7^7 , (N » 1, p > q/2) (48) 

i.e., it is proportional to the system size. Consequently, correlations on a distance k <C N decay like 
A' k/N . This type of behavior is completely different from the one observed in the p < q/2-phase. 
It indicates that the system behaves like the one undergoing sublattice-parallel dynamics |^6[ where 
there is a coexistence phase with non-exponential decay of correlations. 



7 Conclusion 



We have shown that MPS techniques can be applied successfully to an asymmetric exclusion pro- 
cesses on a ring with a defect. In order to apply this method, the model on a ring has to be 
reformulated as a four-state model on a linear chain. This ansatz takes into account that trans- 
lational invariance is broken by the defect while the system is still symmetric under reflections 
with respect to the defect combined with particle-hole exchange. The algebra we derived has very 
simple two-dimensional representations which allow physical quantities like the density profile to 
be computed directly. Another advantage of the matrix formalism is that the model with differ- 
ent dynamical rules can be solved within the same framework (i.e. with the same algebra and 
representations) . 
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The ASEP discussed in this paper is a very simple one where particles hop deterministically in 
the bulk. Our hope is that the present work may show a way how to solve the more general problem 
with probabilistic hopping in the bulk. This generalized model is controlled by three quantities, the 
bulk hopping rate p, the defect hopping rate q, and the particle density p = M/L. In the so-called 
Hamiltonian limit p, q — > 0, p/q = const, it includes the case of random sequential updates which is 
an outstanding problem. 

When introducing a bulk hopping rate p, the algebra (|l8|)-(|l0|) has to be replaced by the bulk 
equations 



[Ai,Bi] 


= 


(i = 1...4) 


(l-p)A i B j 


= BiAj 


€ {(0,2), (1,0), (1,3), (3, 2)} 


AjBi+pAiBj 


= BjAi 


e {(0,2), (1,0), (1,3), (3,2)} 


(1 - p)A i B j + p(l - p)A x B 2 


= BiAj 


G {(0,3), (3,0)} 


(l-p) 2 A x B 2 


= B X A 2 




A 2 B X + pA B 3 + p 2 A 1 B 2 + pA 3 B 


= B 2 A X 





(49) 



and the boundary conditions 

(W\B = (W\A Ao\V) = Bq\V) 

(W\(B 1 + qB 2 ) = (W\A X (l-p)A 1 \V) = B X \V) (50) 

(W\(l-q)B 2 = (W\A 2 (pAx + A 2 )\V) = B 2 \V) 

{W\B 3 = (W\A 3 A 3 \V) = B 3 \V) 

This algebra is even more complicated as the previous one and we were not able to find representa- 



tions or to prove its consistency. However, following the ideas of |19[], one could again assume the 



additional relations given in Eq. ( |22| ) hold. This reduces the algebra to 

pAtAo = A , pA x A 3 = A 3 , pA A 2 = (1 - p)A , pA 3 A 2 = (1 - p)A 3 , 

1 — n (1-p) 2 1 

A A 3 = A 3 A = — — *L-(A 1 +A 2 ), A l A 2 = ±- — ^-A l + — -A 2 , (51) 

u 6 J u p{2-p) y 1 2!) p(2-p) p(2-p) 2 ' y J 

together with the boundary equations 

(W\A 2 = 1—±(W\, A X \V) = -\V). (52) 

Likewise, although this algebra is much simpler, we were not able to prove its consistency and 
the existence nontrivial representations. As usual in ASEPs including a Hamiltonian limit, such 
nontrivial representations are expected to be infinite dimensional. The ASEP with a defect and 
random sequential updates remains as an open problem. 



Acknowledgments 

We would like to thank I. Peschel and in particular A. Honecker for useful hints and interesting dis- 
cussions. S.S. greatfully acknowledges financial support by the Deutsche Forschungsgemeinschaft. 
H.H. would like to thank the Minerva foundation for financial support. 



11 



References 

B. Schmittmann and R. P. K. Zia, Statistical Mechanics of Driven Diffusive Systems, in Phase Transi- 
tions and Critical Phenomena, vol. 17, ed. by. Domb and Lebowitz, Academic Press, U.K. (1995) 

J. Krug and H. Spohn, Kinetic roughening of growing surfaces, in Solids Far from Eguilibrium, ed. C. 
Godreche, Cambridge University Press (1991) 

C. T. Macdonald, J. H. Gibbs, Biopolymers 6 (1968) 1 

M. Schreckenberg, A. Schadschneider, K. Nagel and M. Ito, Phys. Rev. E 51 (1995) 2339 

B. Derrida, E. Domany and D. Mukamel, J. Stat Phys. 69 (1992) 667 
G. Schiitz and E. Domany, J. Stat Phys. 72 (1993) 277 

C. Boldrighini, C. Cosimi, S. Frigio and M. Grasso Nunes, J. Stat. Phys. 55 (1989) 611 
B. Derrida, M. R. Evans, V. Hakim and V. Pasquier, J. Phys. A: Math. Gen. 26 (1993) 1493 
S. Sandow, Phys. Rev. E 50 (1994) 2660 

F. H. L. Essler and V. Rittenberg, J. Phys. A: Math. Gen. 29 (1996) 3375 

L. H. Gwa and H. Spohn, Phys. Rev. A 46 (1992) 844; 
B. Derrida, M. R. Evans and K. Mallick, J. Stat. Phys. 79 (1995) 833; 

G. M. Schiitz, Exact Solution of the Master Equation for the Asymmetric Exclusion Process, to appear 
in J. Stat. Phys. 

V. Hakim and J. P. Nadal, J. Phys. A: Math. Gen. 16 (1983) L213 

A. Kliimper, A. Schadschneider and J. Zittartz, J. Phys. A: Math. Gen. 24 (1991) L955 
and Europhys. Lett. 24 (1993) 293; 

M. Fannes, B. Nachtergaele and R. F. Werner, Europhys. Lett. 10 (1989) 633 
M. R. Evans, D. P. Foster, C. Godreche and D. Mukamel, Phys. Rev. Lett. 74 (1995) 208 

B. Derrida, S. A. Janowsky, J. L. Lebowitz and E. R. Speer, Europhys. Lett. 22 (1993) 651 and 
J. Stat. Phys. 73 (1993) 831 

R. B. Stinchcombe and G. M. Schiitz, Phys. Rev. Lett. 75 (1995) 140 

H. Hinrichsen, S. Sandow and I. Peschel, J. Phys. A: Math. Gen. 29 (1996) 2643 
H. Hinrichsen, J. Phys. A: Math. Gen. 29 (1996) 3659 

N. Rajewsky, A. Schadschneider and M. Schreckenberg, J. Phys. A: Math. Gen. 29 (1996) L305 

A. Honecker and I. Peschel, Matrix- Product States for a One- Dimensional Lattice Gas with Parallel 
Dynamics, preprint cond-mat / 96060*5^ , submitted to J. Stat. Phys. 



K. Krcbs and S. Sandow, Mat rix- Product Eigenst ates for One- Dimensional Stachastic Models and 



Quantum Spin Chains, preprint cond- mat / 9 6 1 00 2 £ 
K. Mallick, J. Phys. A: Math. Gen. 29 (1996) 5373 
M. R. Evans, Europhys. Lett. 36 (1996) 13; 

J. Krug and P. A. Ferrari, J. Phys. A: Math. Gen. 29 (1996) L465 

S. A. Janowsky and J. L. Lebowitz, Phys. Rev. A 45 (1992) 618 and 
J. Stat. Phys. 77 (1994) 35 

G. Schiitz, Phys. Rev. E 47 (1993) 4265 

G. Schiitz, J. Stat. Phys. 71 (1993) 471 

A. Honecker, private communication 



12 



